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Doppler tracking by the DSN is the primary radio metric data type used by 
navigation to determine the orbit of a spacecraft. The accuracy normally attributed 
to orbits determined exclusively with Doppler data is about 0.5 microradians in 
geocentric angle. Recently, the Doppler measurement system has evolved to a high 
degree of precision primarily because of tracking at X-band frequencies (7.2 to 
8.5 GHz). However , the orbit determination system has not been able to fully 
utilize this improved measurement accuracy because of calibration errors associated 
with trzmsmission media, the location of tracking stations on the Earth’s surface, 
the orientation of the Earth as an observing platform , and timekeeping. With 
the introduction of Global Positioning System (GPS) data , it may be possible to 
remove a significant error associated with the troposphere. In this article, the effect 
of various calibration errors associated with transmission media, Earth platform 
parameters , and clocks are examined. With the introduction of GPS calibrations , 
it is predicted that a Doppler tracking accuracy of 0.05 microradians is achievable. 


I. Introduction 

The Doppler data type provides a measure of the line- 
of-sight range rate between a tracking station and a space- 
craft. This functional definition is useful for analyzing or- 
bit determination errors that are spacecraft or trajectory 
dependent, but is of little use for analyzing error sources 
close to the actual measurement such as media or instru- 
mentation errors. The actual measurement is a count of 
a signal derived from the signal received from the space- 
craft and a frequency standard maintained at the tracking 
station that is related to the transmitted signal. Thus, a 
precision model of the Doppler observable would include 


a model of the signal path and hardware elements as well 
as spacecraft dynamics. In practice, the instrumentation 
errors are small compared to media, station location, and 
spacecraft dynamics errors and, therefore, a simple func- 
tional model of the hardware should suffice. 

A model of the Doppler observable is developed below 
that idealizes some of the instrumentation error sources yet 
precisely models the external environment. This model is 
sufficiently precise for computation of the observable and is 
essentially the model contained in JPL’s Orbit Determina- 
tion Program (ODP) [1]. Of particular interest are models 


7 



that are external to the DSN tracking station hardware, 
yet pertain directly to the signal path. Media effects and 
the effect of general relativity on the station clocks are ex- 
amples. Other models, such as station locations and polar 
motion, though not directly part of the Doppler measure- 
ment system, are treated as measurement calibrations for 
the purposes of analysis and orbit determination. 


II. Doppler Measurement Model 

The Doppler measurement is simply an electronic count 
of the number of cycles from a frequency standard ( N c ) mi- 
nus the number of cycles of the spacecraft signal received 
by the ground station (N r ), scaled by the count time in- 
terval ( AT C ). Thus 


(N c - Nr) + n 
AT C 


( 1 ) 


where n is the measurement noise, which is typically about 
1/10 of a cycle. The received frequency and standard fre- 
quency need not be counted individually and differenced, 
but are added together electronically and the beat fre- 
quency counted. This is a detail that, is dependent on the 
hardware implementation. The numerical value of Z m is 
the number that is recorded on the tracking data file and 
is sent to the ODP for orbit determination; the units of 
Zm are hertz. 

In the ODP, one needs to obtain a computed value for 
Zm as a function of parameters that are available to the 
ODP. This function can be derived from the equations 
of motion and the physical model of the system or can 
be worked backwards from Eq. (1) for the measurement. 
Working backwards from the measurement equation is cho- 
sen because the starting point is an equation that repre- 
sents the real physical measurement. 

The frequency standard is obtained by scaling the refer- 
ence oscillator frequency ( f q ) to equal the transmitted fre- 
quency times the spacecraft turnaround ratio which would 
nominally be the received frequency if there were no space- 
craft Doppler shift or additional delay. 

N c = C 3 ft AT C (2) 

where for S-band (2.1 to 2.3 GHz) Doppler, 
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ft = 96 f q 


AT C — T 3e - T 3$ 


The count time (A T c ) is defined as the difference between 
the reception time at the start of the count-time inter- 
val ) and the reception time at the end of the inter- 
val (T 3e ). For a schematic representation of these times, 
see Fig. 1. In Eq. (2), all the parameters are constant or 
arbitrarily specified, including the reception times. The 
real information content of the measurement is contained 
within the count N r . Thus, to obtain a complete equa- 
tion for the computed measurement, one needs an equa- 
tion for N r . It is tempting to differentiate and work in 
the frequency domain; however, the actual observable is 
accumulated phase change, and this is a sufficient reason 
to keep this representation. Therefore, it is necessary to 
formulate the data type in terms of phase thus bypass- 
ing an explicit equation for the received frequency. The 
equation that relates the measurement to the observable 
parameters is simply 


Nr = C 3 Nt 


(3) 


where 


N t = ft (Tu-T u ) 


Equation (3) for N r states that the number of cycles 
counted at the receiver is equal to the number of cycles 
transmitted (N t ) times the spacecraft turnaround ratio. 
This equation is true because the transmitted cycles are 
effectively the same as the received cycles. Thus, the in- 
formation content of the measurement is now contained in 
the transmit times T\ e and T\ s . Since both of these times 
are unknown, some additional equations are needed to tie 
into the observable quantities. At this point in the devel- 
opment, the following equation for the computed measure- 
ment is introduced: 


Z c - (T 3e - T 3s - T\ e + Ti s ) 


Czft 
A T c 


(4) 


Equations are needed for the times in Eq. (4) and these 
will be developed as functions of ephemeris time t. The 
equation for the atomic clock at the station is 
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T = t + F(t,x t y) 


(5) 



p dtdt 


( 10 ) 


where x refers to state variables and y refers to various 
constant parameters. 

The station time ( T ) is equal to the ephemeris time t 
modified by a small correction because of relativity and 
any other parameter that may affect the running of the 
clock. The calibration function (F) is a function of time, 
the solar system gravitational potential, and other con- 
stant parameters. The relevant times relating to the 
Doppler measurement are 

Ti s = ti , + F(ti,,x,y) 

T le = t u + F(t u ,x,y) 

T 3 , = < 3 , + F(t 3s ,x,y) 



are obtained by integrating the equations of motion for 
the spacecraft. These equations are referred to as the light 
time equations and are solved iteratively for the arguments 
of integration. The media delay is included in the measure- 
ment equation by evaluating the calibration function (G) 
at the appropriate times: 


T$ e = ts € -f F(t's ei x,y) 


At” 1 = G(t,x,y) 


(13) 


Making the above substitutions, the equation for the com- 
puted observable becomes anc ^ 


Zc = (<3 e - <3 S - <le + Gj ) ~^r 


A<™ = G(t u ,x,y) 


+ [F(t 3 e ,x,y) - F(t 3s ,x,y) - F(t le ,x,y) 


A tT e =G(t lt ,x,y) 


+F(t i,,x,y)] 


C-ift 

AT, 


( 6 ) 


Since the speed of light is constant in any reference frame, 
one obtains, by integrating along the light path, 


< 3e - <le = +/>23e + A iZ + A it (7) 


At™=G(t 3s ,x,y) 


= G(t 3e ,x, y ) 


The final equation [Eq. (14)] for Z c (the computed mea- 
surement) includes the observable equations as well as 
clock and media calibration functions: 


t 3s - t ls = p Jh±P2h + A<,’’ + At™ (8) 


Pl2 e + P23 e - P 12 5 ~ P23, C^ft 
C ATr 


where the p terms represent the integrated distance along 
the light path and the t m terms represent the additional 
delay caused by transmission media. The distances along 
the light path, 


+ [T(< 3e ,x,j/) - F(t 3s ,x,y) - F(t u ,x,y) 


+Htu ,^.2/)] + [<?(<3e , x, y) - G(t 3s ,x,y) 


P\ 2, 


SC " 


dtdt 


( 9 ) 


+G(<i e ,x,y) - G(ti, , x, y)} 


Gift 

AT, 


(14) 
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III. Analytic Doppler Error Model 

Equation (14) for the computed measurement, while 
useful for actual orbit determination operations, is some- 
what cumbersome for error analysis. An error analysis 
could be performed by combining the sensitivities of each 
term in Eq. (14). An alternative method would be to de- 
velop an approximate analytic representation of the com- 
puted measurement which will be pursued below. The 
range from the tracking station to the spacecraft may be 
described by the following function obtained by integrat- 
ing the equations of motion: 

p=R(t,x,y) (15) 

The velocity of the spacecraft projected along the Earth 
line of sight is obtained by differentiating the range and is 
given by 


In the limit as At approaches zero, this is the definition 
of the second derivative. Without loss of generality, one 
can evaluate the second derivative using two separate ACs 
(one for the first derivative and the other for the second 
derivative) 

<9 2 Z _ F\t + A< 2 + Atj)-F(< + A< 2 ) 
dt* A<jA<2 

_ F(t + AU) - F(t) 

At l At 2 ( 

Solving for the values of the function at the indicated times 
and rearranging the terms gives 

F(t + Ati + A^2 ) ~ F(t + A^) — F(i T A^i) 


9 ~ 


dR 

dt 


(16) 


+ F (0 = Ct a *i A< 2 (23) 


which may be approximated by 

. _ R(t e ,x,y) - R(t, , x,y) 

PK AT c 

From the geometry shown on Fig. 1, then 


(17) 


If one makes Ati equal to the count time (AT C ) and 
A ^2 equal to the round-trip light time (AJ r m), the above 
calibrations can be approximated by the second derivative 
of the calibration function, provided this function has fre- 
quency components with periods that are long compared 
to the delta times. Thus, for the clock calibration function 
evaluations, 


R(te,x,y)-R(t,,x,y) = - (pn t +P23 e -pn, ~Pi3,) (18) 


Substituting Eq. (18) into the equation for the measure- 
ment, the following expression is obtained for the geomet- 
ric part of the computed measurement: 


2 C 3 f t dR 
c dt 


(19) 


Consider the following generic function for the clock 
and media calibration functions: 


Z = F(t) 


(20) 


F(t 3e )-F(t 3s )-F(t u ) + F(t ls ) 


d 2 F(t) 
dt 


2 AT c At r «t 

(24) 


For the media calibration function, one may perform a 
similar approximation. However, because of the reversal 
of signs in Eq. (14), the dominant term is given by the first 
derivative and then 


G(t 3e ) — G(t 3s ) + G(t\ e ) - G(ti s ) & 2 — l. le ^ A T c 


dt 


(25) 


For small A/, the second derivative may be approximated 

by 


Substituting into the exact equation for the measurement 
gives the following approximation: 


d 2 Z F(t + 2A t) - 2 F(t + At) 4- F(t) 


dt 2 


At 2 


( 21 ) 


C 3 ft 


^ + ^i 41rl „ +2 ^ 

c dt 1 dt 


(26) 
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However, if the round-trip light time is too long relative to 
the count time, as would be the case for the outer planets, 
then 


Czft 


2 p dF(t 3 ,x,y) dF(t u x,y) 

c + dt dt 


dG(t 3 ,x,y) dG{t\,x,y) 
H r 


dt 


di 


(27) 


<J Z « I Czft A trtit r, <jj] sin(w e < + A, + <f>) \ <r ae (30) 
where, at X-band and for a typical round-trip light time, 

T r tit = 1, 512 sec 
880 


C 3 = 


749 


The evaluation of the partial derivatives is greatly fa- 
cilitated by breaking up the calibration functions into in- 
dividual terms and approximating with simple functions 
(preferably sines, cosines, or exponentials). The sensitiv- 
ity of the measurement to various ;/ parameters that are 
used to describe the calibration functions is given by 


f t = 32/, + 0.65 x 10 10 Hz 


/, = 20.98 x 10 6 Hz 


<r ae = 3.17679 x 1(T 10 sec/km 


AZ„ 


C 3 J t 


2 + 2 d' 2 G{t,x, y) 


dydt 


dydt 


d 3 F(t,a,p) 

dydt 2 


A i r 


tit 


A y 


(28) 


This equation gives a peak amplitude of 70 mHz for a 100- 
percent error in the amplitude (a e ). Thus, to obtain a 
1-mHz accuracy in the measurement, the calibration error 
must be less than 2 percent. 


The y parameters in Eq. (28) can be formed into a co- 
variance matrix which, when pre- and post-multiplied by 
the measurement sensitivities, gives the measurement vari- 
ance. 

For example, consider the case where the only y param- 
eter is the amplitude of the daily special relativity term 
that is part of the atomic clock calibration function. This 
term is defined by 

F e = ot e r s sin (u € t + X s + 0) (29) 


IV. Data Noise Evaluation 

Approximately two orbits of Magellan data were pro- 
cessed from February 7, 1991, 19:29:02 (ET) to February 
8, 1991, 00:10:33 (ET). The postfit residuals are shown 
on Fig. 2, which involved solving for the spacecraft state, 
a constant nongravitational acceleration, and a fourth- 
degree and -order spherical harmonic gravity field. The 
first 10,000 sec of data shown on Fig. 2 are from Deep 
Space Station (DSS) 15. The gap around 6,000 sec is data 
unavailable near periapsis. Data after 12,000 sec are from 
DSS 45. 


where a e is the subject y parameter, r $ is the station spin 
radius, u e is the Earth rotation rate, A, is the tracking 
station longitude, and </> is a phase angle needed to bring 
UT1 into agreement with ephemeris time. Typical nominal 
values of these parameters are given by 


Recall from Eq. (1) that the measurement noise is 


N c — N r n 

-atT + at c 


(31) 


a c = 3.17679 x 10” 10 sec/km 
r s = 5,204 km 
u) e = 7.292 x 10~ 5 rad/sec 

Taking the required first and second partial derivatives, 
one obtains the following equation for the measurement 
error as a function of time: 


The data noise is approximately 1/10 of a cycle count, and 
is independent of frequency. For a 60-sec count time, the 
Doppler data noise is therefore about 1.66 mHz. This value 
agrees with Fig. 2, which shows noise of this magnitude 
from DSS 15 and appears to be uniformly distributed but 
with some structure. The data from DSS 45 are a bit more 
noisy and appear more Gaussian. 

Doppler data are scaled by the count time to make the 
recorded measurement proportional to range rate. The 
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Doppler measurement sensitivity to line-of-sight velocity 
is given by Eq. (19): 

p (32 ) 

c 

At S-band, typical values for the constants in the above 
equation are: 

240 
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ft = 96 f q 
f q = 22 x 10 6 Hz 
c = 299792.458 km/sec 


Rd = fd(zd, t) 

(36) 

Zd = 

(37) 

Rw - fd{Zu jj7) 

(38) 

z w = 

(39) 

■*0 

< 

II 

(40) 


In taking the required partial derivatives with respect to 
time, one assumes R w and Rd are linear in z w and Zd and 
7 is a function of only time: 

dG t dR d . OR-ddzd t dRu . 0R w dz w /yll . 

■ar = _ a7 1 ' + J^-m + W y+ 0771T (41) 


and solving Eq. (32) for p gives 0.108 mm/sec for Zp = 
1.66 mHz data noise. One can also obtain from Eq. (32) 
the well-known result that for p = 1 mm/sec the measure- 
ment noise is Zp — 15.3 mHz. At X-band, the measure- 
ment noise Z p — 1.66 mHz corrresponds to 0.03 mm/sec 
velocity. 


The expression for the measurement sensitivity becomes 


AZ„ 


2C 3 /« 


(P_Rd . dRd d^zj 

[dyd~( 7 + dz d dydt 


, d 2 R w . dR w d 2 z u 
+ — ^ 7 + 


dydy dz w dydt J 


Ay (42) 


V. Troposphere Calibration 

A radio signal passing through the Earth’s troposphere 
will be delayed depending on the dielectric constant of the 
media and path length. A functional formula for this delay 
has been defined above. 

At 1 = G t (t y x> y) (33) 


The first term in the above equation represents the non- 
linearity of the dry troposphere mapping function, and 
the second term represents the variation in the dry tropo- 
sphere z height due to local weather. The next two terms 
are the same quantities for the wet troposphere. The tro- 
posphere wet and dry mapping functions are tabulated as 
delay as a function of spacecraft elevation angle. Empirical 
formulas for these mapping functions are given by [2,3] 


The sensitivity of the Doppler measurement to media 
is functionally given by Eq. (28) and one may extract the 
part that pertains to the troposphere. 


Rd — 


Zd 


C 


sin 7 -f 


Ad 


Bd + tan 7 


AZ rn 




(34) 


The troposphere delay has been separated into wet and 
dry components that are functions of delay at zenith {z) 
and elevation angle (7): where 




sin 7 4- 


B w + tan 7 


(43) 


(44) 


G t (t,x,y) = R d + Rw 


(35) sin 7 = cos 6 cos A cos <j> + sin A cos <j> + sin </> sin £ (45) 
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A — -f- A j — ct 


(46) The partial derivatives of 7 and z w with respect to time 
are given by 


The dry component of the troposphere (Rd) is a function 
of the delay at zenith (z^), the elevation angle ( 7 ), and 
constants Ad and Bd that are provided to model the bend- 
ing at low-elevation angles. The wet component ( R tu ) is 
similarly defined. The elevation angle ( 7 ) is computed as 
a function of the latitude of the tracking station (</>), the 
declination of the spacecraft ((5), and the local hour angle 
with respect to the spacecraft (A). The local hour angle 
is zero when the spacecraft is at zenith and is a function 
of the Greenwich hour angle (c J e t), the station longitude 
(A 5 ), and the right ascension of the spacecraft (a). 

The dry component of the troposphere is assumed to 
be stable and most of the variability is associated with 
the wet component. The variation in the wet component 
may be modeled as a periodic variation in the z height 
( z w ). The actual variation obtained by tracking the GPS 
satellites in early 1991 is shown on Fig. 3. The hourly 
variation in the wet component of the troposphere appears 
as a random walk that would require a high-order Fourier 
series to represent analytically. For the purpose of error 
analysis, the variation may be modeled as a simple sinusoid 
with amplitude and frequency selected to be representative 
of Fig. 3. 


cos <j)(cos A — cos 6 sin A) • 

A 

(51) 

cos 7 

A = uj e 

(52) 

Zw = Z Wl U> Wl COS (Wu,^) 

(53) 


One is interested in the sensitivity of the Doppler measure- 
ment to the constant and periodic terms of the expression 
for the wet troposphere z height: 


AZ m 


2C 3 f 1 d 2 R w . A 
— d^ 7Ay 


(54) 


Substituting the above partial derivatives into this equa- 
tion yields 


A Z 20 


2 C 3 f t 


cos 7 7 
sin 2 7 


Az< 


t ^ 0 


(55) 


z u> - z WQ + Z WI sin (w Wl t) 


(47) 


For the troposphere, the y parameters that are of inter- 
est are the constant coefficient of the wet z height function 
(z WQ ) and the periodic coefficient (z Wx ). In order to obtain 
these sensitivities, it is necessary to evaluate the partial 
derivatives contained in the above measurement sensitiv- 
ity function [Eq. (59)] that pertain to the wet troposphere 
component. For elevation angles above the horizon by a 
few degrees, one can approximate the mapping function 
by 



Z\d 


c sin 7 


and 


dR w 1 — z w cos 7 

dy c sin 2 7 


(48) 


(49) 


and so 


AZ z 


2C 3 f t 


x 


/ 1 

\ — w Wl COS t - 

f sin 7 1 1 


cos 7 7 

, 9 sin ^ xd . t 
sin 7 1 


A z w ^ 


(56) 


An example of the effect of the troposphere on X-band 
Doppler is the application of GPS troposphere calibration 
data to the Magellan spacecraft. The GPS data were ob- 
tained at a time when the spacecraft was in orbit about 
Venus. It would be preferable to have cruise data, where 
the effect of the troposphere is more pronounced, rather 
than orbiting data where separation of the tropospheric 
effect from the gravity field modelling effects is difficult. 
Recall the two orbits of Magellan Doppler residuals shown 
on Fig. 2. At the time of the first data point, the param- 
eters in the above Doppler sensitivity equations had the 
following values: 


dRxj) _ 1 1 

dz w c sin 7 


(50) 


C 3 = 


880 

749 
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7 = 37.49 deg 


An empirical formula for the effect of the ionosphere on 
the Doppler measurement is given by [4] 


ft = 32 f q 4* 0.65 x 10 10 Hz 
7 = 3.818 x 10 -5 rad/sec 


Gi = —J2kCjXl (59) 

c i= o 


f q = 20.98 x 10 6 Hz 
z WQ — 5 x 10“ 5 km 
Z u , 1 = 5 X 10“ 5 km 
w Wl = 1.454 x 10~ 4 rad/sec 


- 1 <6o) 

where the C/s are coefficients of a polynomial in time ( t ) 
from t a to 4 normalized over the interval of -1 to -hi and 
k is a proportionality factor introduced for the purpose of 
computing error sensitivity. The first partial derivative of 
Gi with respect to t is needed: 


The above analytic Doppler sensitivity equations re- 
veal a Doppler measurement error sensitivity of 0.044 mHz 
per centimeter of constant wet troposphere z height (zu, 0 ) 
error and a periodic Doppler measurement sensitivity of 
0.044 mHz in quadrature with 0.134 mHz per centimeter 
of variable wet troposphere z height amplitude (z^). Fig- 
ure 4 shows the prefit residuals for the periodic wet tro- 
posphere 2 height calibration defined above. A Doppler 
shift of about 0.67 mHz is apparent early on for the data 
from DSS 15; at the end of Fig. 4, the Doppler shift is 
about 2 mHz for the low elevation data from DSS 45. The 
data in Fig. 4 indicate that the error contribution of the 
troposphere, which may amount to several mHz, could be 
reduced or eliminated with GPS calibration data. 


VI. Ionosphere Calibration 

A radio signal passing through the ionosphere experi- 
ences a reduction in group velocity and an equal increase 
in phase velocity that is a function of the frequency and 
the number of charged particles along the signal path. The 
Doppler measurement is dependent on the phase velocity 
and the advance of the signal is functionally defined by 


0Gi = _2fc_£ (gl) 

dt 3 

The sensitivity of the Doppler measurement to the con- 
stant of proportionality k is given by 

AZ k « J2 j C > Xj 1 Ak ( 62 ) 

C Tfc ta j — i 

For the first data point shown on Fig. 2, the following val- 
ues for the constant parameters are needed to compute the 
Doppler sensitivity in addition to some of the parameters 
given above for the troposphere: 

t = -280,772,197 sec 
t a = -280,785,442 sec 
t b = -280,745,040 sec 


At* = Gi(t,x,y) (57) 


Cj =5.8432, -1.2655, 5.8368, -1.0912, 


The sensitivity of the Doppler measurement to media is 
functionally given by Eq. (28) and may be adapted to the 
ionosphere: 


A Z m 


2 Csft 


dydt y 


(58) 


- 0.6015, 0.5551 m 


Assuming a 10-percent error in the ionosphere model 
(Ak — 0.1), the error in the Doppler measurement at- 
tributable to the ionosphere is 0.143 mHz. 
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VII. Station Location Errors 

The Doppler measurement of the distance from a track- 
ing station antenna to a spacecraft enables the orbit deter- 
mination system to estimate spacecraft state, station loca- 
tions, planetary ephemerides, and many other parameters 
that are of interest for navigation. If the primary pur- 
pose of the estimation process is determining the orbit of 
a spacecraft, then the station locations and the orientation 
of the Earth’s crust in inertial space are often regarded as 
part of the measurement system. Station locations, timing 
and polar motion, continental drift, and solid Earth tides, 
to name a few quantities, are calibrated out of the data. 

Regarding the station locations as part of the measure- 
ment system, the following function may be used for error 
analysis and approximates the effect of the tracking station 
when the spacecraft is a great distance from the Earth: 

R s — - r s cos 6 s\n(u> e t + A, ) (63) 

c 

This function can be regarded as a calibration of the data 
and the Doppler measurement sensitivity is functionally 
the same as for the troposphere and ionosphere. The range 
is premultiplied by one over c to give this function the units 
of time consistent with the other calibration functions. 


AZ m 


2C 3 f t 


d 2 R s (t, x, ?/) 
dydt 


Ay 


(64) 


At present, the DSN station locations have been deter- 
mined in the terrestrial (Earth-fixed) frame with an ac- 
curacy of about 0.1 m. Assuming a 0.1-m error in the 
tracking station spin radius (r 5 ), the error in the Doppler 
measurement attributable to station locations is about 
0.4 mHz. 


VIII. Clock Calibration 

According to the theory of relativity, a clock running 
in a frame of reference that is moving with respect to an 
observer’s frame of reference or in a gravitational potential 
field will appear to run slower by an observer that is as- 
sumed stationary and is removed from the potential field. 
Therefore, a distant observer who is stationary with re- 
spect to the solar system will see the atomic clocks at the 
DSN stations running slower than his or her hypothetical 
clock. The observer’s clock records coordinate time, which 
is called post-Newtonian time (PNT) [5,6]. 

The relationship between PNT and the proper time 
measured by an atomic clock is given by the metric. For 
a particle moving in an orbit around the Sun, the metric 
in isotopic Schwarzschild coordinates is given by 

ds 2 — ^1 c 2 dt 2 — v 2 dt 2 -f 0(c~ 2 ) 


where 


The sensitivity of the Doppler measurement to the dis- 
tance of the tracking station from the Earth’s spin axis 

(r s ) is: 



A Z r$ « u> e cos 6 cos (w € t -p A 5 ) A r s (65) 

c 

For the first data point shown on Fig. 2, the following 
values for the constant parameters are needed to compute 
the Doppler sensitivity in addition t.o the parameters given 
above for the troposphere and ionosphere: 


r s = 5, 204 km 


Solving for proper time (ds 2 = c 2 dr 2 ) yields 



which can be further approximated by 

dr _ ^ It ; 2 /x e 

dt c 2 r 2 c 2 c 2 r e 


( 66 ) 


(67) 


( 68 ) 


cj c = 7.292 x 10 5 rad/sec 


6 = 8.39 deg 


where the Earth’s gravitational potential is separated from 
the Sun. The atomic clock time (r) is obtained as a func- 
tion of t by integrating the metric in conjunction with the 
equations of motion: 
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which can be evaluated directly by numerical integration, 
< 69 > giving 


where 


r _ ill. _l i 

c 2 r + 2 c 2 c 2 r e 


r=i+ C~ L (t-^ + P-H‘' <(77) 

( 7 °) 

This equation will be referred to as 


The function 1/ can be separated into a constant term (Lo) } 
secular terms that grow with time (L s ), and periodic terms 
(i P ). Thus 


Moyer [1] provides an approximate analytic formula for 
the periodic terms, attributing much of the original work 
to Brooks Thomas: 


L — Lo + L s -f L p (71) 

The constant term ( L 0 ) is obtained by averaging L over 
all time and can be represented by 


T « t - L 0 (t - to) - — (rj ■ rj) - -^(fj • r*) 


— fr c • r e> ) - — ('r c • rH — (r- • r-) 

c* (e r «> c 2^» r <>) c 2 (/) . +(it )'j r J> 



— + o ' 0 

K r 0 Z 


) 


+ 


c 2 r e 


(72) 


f*sa 


C 2 (n,a +/i s ) 


(r a ■ r a ) 

4 ia/ 


(78) 


where ro and i>o are constants that give the correct value of 
Lo in the above equation. For the Earth’s orbit about the 
Sun, r 0 is approximately the semimajor axis of the orbit 
and i s approximately the mean orbital velocity. Since 
the orbit is nearly an ellipse, 

- V 2 (73) 

a r 

and for r = a, 


In the notation used above, the position of the body identi- 
fied by the subscript is with respect to the body identified 
by the superscript, where c is the solar system barycen- 
ter, s is the Sun, b is the Earth-Moon bary center, e is the 
Earth, j is Jupiter, and sa is Saturn. This equation will 
be referred to as the vector function and is designated by 

F v . 

Another formula is provided in part 2 of [1] that de- 
scribes the periodic terms as functions of sines and cosines. 
The dominant terms are given by 


Lo 


3 ji s 

2c 2 a 


(74) r — Lo (t — fo) - 1-658 x 10 3 sin E 


The secular terms L , are assumed to be zero because of 
conservation of energy and momentum. This leaves the 
periodic terms and these are given by 


- 1.548 x 10 -6 sin D 

- 3.17679 x 10 -lo r a sin (UT1 + A) - ... (79) 


-? O-s+P-H 


(75) 


Recall that the station clock calibration function is defined 

by 


r = t + F(t,x,y) (76) 


where E y D> and UTl are the angles describing the Earth’s 
orbit about the Sun, the Moon’s orbit about the Earth, 
and the rotation of the Earth about its spin axis, respec- 
tively, and A is the ground longitude. If one replaces the 
sines and cosines in the complete expression for r by a 
power series and truncates, there is a polynomial in t . The 
above equation will be referred to as the polynomial func- 
tion and is designated by F p . 
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Since there are three functions that purportedly give 
the same result, it may be instructive to tabulate the dif- 
ferences as a function of time. One can then apply the 
error sensitivity given by Eq. (28) to determine the error 
in the Doppler measurement. This procedure will yield 
the error in the differences between the functions but will 
not reveal which function is correct. Indeed, the correct 
function may be some unknown fourth function. 

As a basis for this comparison, the function (Fi) ob- 
tained by numerical integration is selected. Figure 5 shows 
a plot of UTC minus ephemeris time (ET) as a function 
of ET. The ET is defined as PNT time with the constant 
rate term Lq removed. A detailed discussion of this sub- 
ject is provided by Hellings [6]. Dropping the Lq term 
from the definition of ET results in the length of the ET 
second being shortened. Therefore, in order to compen- 
sate, the masses of all the bodies in the solar system are 
scaled. UTC is simply an atomic clock time r adjusted by 
a constant that includes accumulated leap seconds. The 
1991 value of this constant is —58.184 sec. 

In order to perform the integration to obtain F{, a 
value for Lq and the constant of integration must be ob- 
tained. An initial guess for Lq is obtained from Eq. (41) 
and the constant term is obtained from the vector func- 
tion F v [Eq. (46)]. The equations of motion are integrated 
for 1 year and a mean value of Lq is obtained which con- 
strains the function F{ to equal the vector function F v at 
the end points of the interval. The linear term Lq obtained 
in this manner has a value of 1.55035 x 10~ 8 , which com- 
pares quite favorably with the linear term, given in [2], 
of 1.55052 x 10 -8 . The difference may be attributed to 
the averaging interval or the effect of Jupiter and Saturn 
on the solar system barycenter. Next, the integrated vec- 
tor and polynomial clock functions over the same 1-year 
time interval are evaluated, and the differences F v — F( 
and F p — F{ are plotted on Fig. 6. Inspection of Fig. 6 
reveals an annual error term with an amplitude of about 
1.5 x 10 -6 sec in both the vector and polynomial func- 
tions, a biweekly term of about 1.0 x 10 -7 sec in only the 
polynomial function, and a daily term that is less than 
1.0 x 10 -12 sec. 

The error functions shown on Fig. 6 are quite small 
compared to the actual clock calibration function shown 
on Fig. 5. As previously noted, the error functions arise 
from differences in the methods for computing the func- 
tion and as a result provide very little insight into the 
error source. The error source may be high-order terms 
that have been neglected, approximations that have been 
made, or simple programming errors in implementation of 
the comparisons. However, the biweekly oscillation in the 


polynomial function is highly suspect because this term 
has been omitted from the expansion. Without judging 
the quality of the approximations, one proceeds to com- 
pute the effect of these error terms on the Doppler mea- 
surement. 

An error function can be defined that includes all of the 
suspect terms and this is given by 

Fe = t € + L e (t - t 0 ) + A d sin(u^f) 

-F A bw sin(2u> bw t) + A y sin(u ; y t) (80) 

where A dy A bxi} , and A y are the amplitudes of the daily, 
biweekly, and annual terms identified above, respectively. 
The clock offset error is r e and the clock rate error is L e . 
Recall from Eq. (27) that the sensitivity of the Doppler 
measurement to errors in the amplitude (A) is given by 

AZ m » C 3 ft d3F * { -'A V) A< r(/t AA 
oAot 

and the resulting error is given by 

^2 » \C a f t A trtit w 2 sin(wf)| cr A (81) 

Observe that the constant and linear terms do not con- 
tribute to the measurement error. Numerical values for 
daily, biweekly, and annual terms are given in Table 1. 
The Doppler measurement errors associated with the daily 
and annual terms are negligible. The biweekly error of 
0.0345 mHz, though small for Doppler measurements, may 
present a problem for VLBI measurements. 

IX. Summary 

A survey of the X-band Doppler measurement system 
has supported the generally accepted notion that tropo- 
sphere calibrations are the dominant error source. Using 
only a seasonal troposphere model, calibration errors in 
the troposphere contribute 1 or 2 mHz to the measurement 
error. With the GPS calibration of the troposphere, one 
may be able to obtain a factor of five reduction in the tro- 
posphere error. With station location errors at about the 
10-cm level, the inherent Doppler measurement accuracy 
of 1 to 2 mHz may be fully utilized. A 1-mHz Doppler mea- 
surement error would translate into about a 0.02 mm/sec 
error in the line-of-sight velocity and, according to recent 
studies [7], perhaps 50 nr ad in angular measurement accu- 
racy. 
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Table 1. Doppler error sensitivities to clock error terms. 


Error term 

Amplitude, 

sec 

Frequency, 
rad /sec 

Doppler error, 
Hz 

Daily 

1.0 x io- 12 

7.292 X 10“ 5 

6.77 x 10“ 8 

Biweekly 

1.0 x 10“ 7 

5.209 X 10" 6 

3.45 X \0~ b 

Annual 

1.5 X 10“ 6 

1.991 X 10“ 7 

7.57 x 10 -7 
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